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Abstract. In the manganites REi—x^ExMnOs {RE and AE being rare-earth and alkaline- 
earth elements, respectively) the random distribution of RE'^^ and AE'^^ induces random, but 
correlated shifts of site energies of charge carriers in the Mn sites. We consider a realistic model of 
this diagonal disorder, in addition to the double-exchange hopping disorder, and investigate the 
metal-insulator transition as a function of temperature, across the paramagnetic-ferromagnetic 
line, and as a function of doping x. Contrary to previous results, we find that values of 
parameters, estimated from the electronic structure of the manganites, are not incompatible with 
the possibility of a disorder induced metal to insulator transition accompanying the ferromagnetic 
to paramagnetic transition at intermediate doping {x ~ 0.2 — 0.4). These findings indicate clearly 
that substitutional disorder has to be considered as an important effect when addressing the 
colossal magnetoresistance properties of manganites. 
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1. Introduction 

The discovery of colossal magnetoresistance (CMR) in the manganites i?i?i_a;A£'a;Mn03 (where 
RE and AE are trivalent rare-earth and divalent alkaline-earth ions, respectively) has attracted 
much interest to these perovskites On the one hand, understanding the origin of the CMR 

effect from a fundamental point of view is expected to give some insight into the complex behaviour 
seen in other strongly correlated systems, as high-temperature superconductors. On the other, 
such a colossal response to an external perturbation still makes these Mn oxides very appealing 
from the point of view of applications ^]. This CMR effect, particularly "colossal" in the so-called 
intermediate-bandwidth manganites for doping x ^ 0.2 — 0.4 [4], is intrinsically related with the 
presence of metallic behaviour below the Curie temperature (Tc) and insulating behaviour above 
it. Understanding the nature of this metal-insulator transition (MIT) accompanying the magnetic 
transition is thus a key point in the comprehension of the CMR effect. 

Qualitatively, the correlation between transport and magnetic properties is well understood via 
Zener's double-exchange (DE) mechanism [3-[9|]: the spin of itinerant Cg and local t2g Mn d electrons 
are aligned by Hund's rule; to lower the kinetic energy ferromagnetism is favoured, and at low 
temperature (T < Tc) a ferromagnetic half- metal is realised [7, 'if . Consequently, early proposals 
for the MIT focused on the random nature of hopping in the paramagnetic phase (T > Tc) [10-13J. 
Quantitative analysis based on the transfer matrix method showed concomitant ferro-paramagnetic 
and MITs for x ~ 0.2 — 0.4 when diagonal disorder is added to the DE model [3 llEl- The resultant 
CMR effect |il(i.lT7|. however, requires an amount of diagonal disorder that seems incompatible with 



potential fluctuations originated by the random distribution of RE^^ and AE^^ ions p, ll8| . This 
apparent failure of Anderson localisation theories favoured models based on polaronic formation 



appare 

ii-ii 



2l|, owing to the electron-phonon coupling due to Jahn- Teller effect in manganites [5]. It has 
been argued, however, that manganites fall into an intermediate electron-phonon coupling regime 
where small-polarons - key ingredients for a MIT driven by electron-phonon interaction - are hardly 
formed in the paramagnetic phase [1, [l^l . More recent theories explain the CMR as an effect of 
competing orders: in brief, a ferromagnetic metal competes with a charge-order insulator phase, 
producing, for T > Tg , an inhomogeneous state highly sensitive to external perturbations where 
CMR is observed 23|-|29j. The underlying model producing such a phase competition scenario 



takes into account, on the same footing, the DE mechanism and the electron-lattice coupling and, 
surprisingly, intrinsic disorder in manganites. Adding quenched disorder has been shown to make 
the inhomogeneous state even more sensitive to external perturbations, enhancing the CMR effect, 
and avoiding fine tuning of model parameters (30| . However, the need to include disorder on the 
same level as the DE mechanism and coupling to Jahn- Teller phonons strongly contradicts the 
observation that disorder in manganites is sufficiently weak for a virtual crystal approximation to 
be reasonable 0, ■ 

In this paper we consider a realistic model for diagonal (substitutional) disorder in manganites, 
in addition to the DE hopping disorder, and investigate the MIT as a function of temperature, across 
the paramagnetic-ferromagnetic line, and as a function of doping x. Contrary to previous results 



we find that values of parameters, estimated from the electronic structure of manganites |18ll3l|. are 
not incompatible with the possibility of a disorder induced MIT accompanying the ferromagnetic to 
paramagnetic transition at intermediate doping {x ^ 0.2 — 0.4). Therefore, substitutional disorder 
has to be considered at least on the same foot as the coupling to the lattice when addressing the 
CMR properties of manganites. These findings give support to theories where disorder is a key 
ingredient (32| . as the phase-competition scenario for CMR, where, as mentioned above, diagonal 
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disorder plays a crucial role [23|-|30|. 
2. Model 

In order to model substitutional disorder in manganites we note that, for each carrier introduced 
in the system, there is a RE^'^ — )■ AE'^'^ substitution. The corresponding change in the Coulomb 
field shifts the site energy of an electron in a manganese site at a distance R by 

where e is the relative dielectric constant of the material. We can take this effect into account by 
including a random site energy term in the DE Hamiltonian, 

^ = - + ^■^■) + E (2) 

{ij) i 

The first term on the right hand side in ^ is the usual infinite Hund coupling DE Hamiltonian, 
where the hopping of Cg electrons between nearest-neighbours Mn sites depends on the background 
configuration of classical t2g core spins = 5 (sin 6i cos <f)i , sin 6i sin 0; , cos 6i ) , with 

t(Si, Sj) = t[cos{e,/2) cos(6'j/2) + e-'^"^'-^^) sin(6ii/2) sin(6ij/2)]. (3) 

The second term on the right hand side of ([2]) stands for the diagonal site disorder. It has 
been modelled in previous works with a uniform probability distribution for —1^/2 < Ci < W/2 
(Anderson disorder) (l^ - [l6| . Through the analysis of the mobility edge trajectory in the energy 
vs disorder (W) plane, obtained using the transfer matrix method [33-35], it has been found that 
a MIT occurs when the system crosses the ferro-paramagnetic transition line, for 0.2 < x < 0.5 
provided the diagonal disorder is strong enough, 12t < W < 16.5t. The plausibility of such large 
value of the disorder parameter, however, has been questioned Q mainly on the basis of density 



functional results obtained by Pickett and Singh [18|, |3l| . In the following we show that a thorough 
analysis of the results of and [3l| are, in fact, not incompatible with disorder in the range 
12t<W < IG.bt. 

2.1. Substitutional disorder strength 

Pickett and Singh [H, 31| looked at the a: = 1/3 concentration, and performed LDA calculations 



of band structure for a periodic structure of La2CaMn309 with a tetragonal unit cell containing 
a La-Ca-La set of planes. There are two inequivalent Mn sites in this structure, one with eight 
La^^ and the other with four Ca^"^ and four La'^"'' nearest- neighbours (NNs). The local density of 
states at the Mn sites showed a difference of AeMn ~ 0.5 eV between the band edges for these two 
types of sites, which was interpreted as arising from the different charges Ca^^ and La'^^. Thus, 
AeMn = 4Vl with Vi given by H]) with R = a%/3/2 (the first NN La - Mn distance), 

2e2 

where a « 3.9 A is the Mn — Mn distance {sl]. From the calculated value of AeMn a dielectric 
constant e « 34 is obtained, or equivalently Vi ~ 0.13 eV ~ 0.6t [t « 0.2 eV was used Ql)- 

Such a dielectric constant, however, is quite unlikely. We should note that equation (j4]) is 
actually a microscopic description, where R — a-\/3/2 « 3.4 A. Neglecting metallic screening, we 
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should get a relative permittivity reflecting the polarizability of the Mn d—Op complex, as pointed 
out in |l8| . Infrared reflectivity measurements on Lao.erCao.asMnOs give a high-frequency dielectric 
constant £00 ~ 7.5 at 78 K [37|, which, though being only a lower bound, casts serious doubts on 
e ~ 34. On the other hand, note that the result AeMn — 4Vi is a special case where only first NNs 
(La/Ca sites) contribute to the local potential. A more realistic situation should account for next 
NNs contributions. 

It is easy to generalise the first NN result AeMn = 4Vi in order to account for the Coulomb 
contribution of the the ith shell, Vi = e^/ {AireoeiRi), where £i is the dielectric constant for the given 
shell. In particular, taking into account second and third NNs, we get AeMn = 4Vi — I2V3, where V2 
is absent because the two inequivalent Mn sites have the same second NN environment. The value 
AeMn ~ 0.5 eV found by Pickett and Singh fisj is reproduced with Si ~ 10 and £3 ~ 17, where we 
used i?3 = a\/T9/2 « 8.5 A. Following [18j, we will keep only first and second NN contributions, 
with El — 62 ~ 10, and R2 — avTT/2 f« 6.5 A. The resulting random site energies may be written 
as 



e^ = Vlikl+k2V3/n), (5) 

where lij is the number of AE"^^ ions in the jth shell of Mn site i for a given RE^^ /AE'^'^ 
configuration. Inserting Ex « 10 in equation (U) we get Vi ~ 0.43 eV sa 2 At. The probability 
distribution for this random, but correlated site energy model may be written as 



ii2=o ^ ' hi=a 



(l-x) 



32~lii~li2 I 



(6) 



In figure [T] we show the resulting coarse grained distribution as a function of e^ — (e^) for x = 1/3, 
where the average site energy is (e^) — xVi{8 + 24^3/11); the inset shows the true discrete Mn- 
site energy probability, as given by the weight of the delta functions in ([6]). The distribution is 
approximately Gaussian with a root mean square (RMS) deviation 



2) - (e,)2 ^ AM, (7) 

as obtained by fitting with a Gaussian distribution (dashed line in figure [1]). A rectangular 
distribution with the same RMS deviation has W = a^/T2 w 15. 9t; well in the range required 



for a MIT at Tc, 12t < W < 16. 5i [15|, |17|. The effect on Tq of a discrete random site-energy 
distribution similar to equation ([6]), but restricted to the NN shell, has been considered in (38| . 

The analysis which lead us to figure [1] and equation ([7]) was previously carried out in [3| , 
but with a very different conclusion about the strength of disorder in OMR manganites. The 
main difference with respect to the present analysis is the value of the dielectric constant entering 
equation ^ through V^. Even though it is stated in [3| that first and second shells of RE^^ / AE'^'^ 
sites are taken into account with ei — 62 ~ 10, we can only reproduce figure 1 of [3| (the analogous 
of our figure [1]) if equation ([6]) is used with si — 62 ^ 34 (an unlikely dielectric constant, as 
discussed above). As a consequence, the associated distribution was found to have a RMS of 

f An additional difference comes from the misplaced second shell of RE'^'^ / AE'^'^ sites in llSll . The distances to 
the first and to the second shells differ by 48%, and not by the referred 13% in |l8l |. This error, however, does not 
change considerably the results: while in the case of 48% difference the second shell has 24 sites, in the case of 13% 
difference it has only 6 sites. 




Figure 1. Full line shows the probability distribution of Mn-site energies due to random 
placement of RE^'^ (2/3 probability) and AE'^'^ (1/3 probability) on first and second neighbour 
sites, obtained by substituting (5— functions in ^ by Lorentzians with half width t at half 
maximum. A dielectric constant of e 10 was used for both shells. The dashed line is the 
fit to the full line with a Gaussian distribution. The inset shows the true discrete Mn-site energy 
probability. 



cr « 1.3< « 0.26 eV[|. This means that a rectangular distribution with the same RMS deviation has 
W = aVT2 « AM; well below the values required for a MIT at Tc, Ut<W< 16.5t 

3. Results and discussion 

From the above considerations we may conclude that a realistic parametrisation of diagonal disorder 
is mandatory for a precise estimate of the effect of substitutional disorder in manganites. We have 
performed a transfer matrix calculation (33l - [35| using the model given in equation ([5]), with site 
energies calculated from a random distribution of the dopant ions AE'^^. First and second shells 
of RE^^ / AE^^ -sites were taken into account assuming equal dielectric constant, which results in 
random on-site energies given by ([5]), and a site energy distribution given by In this case 
the site disorder is parametrised by x, which determines the fraction of AE'^'^ ions in the system 
and thus the variables lij in [SI and by the parameter Vi given in equation |31 or equivalently the 
dielectric constant e. The doping level x also determines the Fermi energy in the system, which 
was calculated by integrating over the disorder averaged density of states obtained for clusters of 
64 X 64 X 64 sites using the recursive Green's function method [39]. 

The transfer matrix technique gives reliable information about the extended or localised nature 
of the eigenstates [s^ 113 • For a quasi-one-dimensional system with length L and cross-section 
M X M, where L ^ M (in units of the lattice constant), the method provides the localisation 
length Am of the finite system from the smallest Lyapunov exponent of the respective transfer 
matrix product. The scaling behaviour of the normalised finite-size localisation length Aaz/M then 

§ In list the full width of the distribution at half maximum was found to be &t ^ 0.6 eV. If we assume the distribution 
to be Gaussian, the full width at half maximum is related with the root mean square a as Se = a2 V21n2. 




Figure 2. Critical value of Vi vs x in the paramagnetic (circles) and ferromagnetic phases 
(triangles). A system with (a;, Vi) in region I (III) is metallic (insulator) in the paramagnetic 
and ferromagnetic phases; in region II the mobility edge crosses the Fermi level in the para- 
ferromagnetic transition. The inset shows the typical behaviour of the normalised finite-size 
localisation length Aj\//A/ vs Vi/t across the MIT for the particular case o{ x = 0.30 in the 
paramagnetic phase. 



determines the nature of the eigenstates for a given x and Vi at the Fermi energy. Extended 
(locahsed) states show increasing (decreasing) Xm/M as M increases. This behaviour can be 
appreciated in the inset of figure [5] for x — 0.30 in the paramagnetic phase, where the direction of 
each spin is chosen randomly from a uniform distribution on a sphere. We have used system bars 
{L X M X M) with a longitudinal length L such that the relative error in Aj\/ is < 1% (typically 
L ^ 10^), and M — 6,8, 10, 12. At criticality Xm/M is independent of M, signalling the Anderson 
transition, and providing the critical parameter values [41]. 

The main result of this work is presented in figure [51 For each concentration the critical 
values of Vi, at which the mobility edge and the Fermi level coincide, were calculated in the 
paramagnetic (circles) and ferromagnetic phases (triangles). A value of Vi between these two 
(i.e.. in region II) implies a crossing of the Fermi level and the mobility edge when the system 
orders magnetically. A value of Vi w 3t is sufficient to give rise to an Anderson MI transition for 
concentrations x ^ 0.2 — 0.5. While it is still higher than the estimate based on [18] and [31i] , 
Vi « 2.1t, it is sufficiently close to cast some doubt on a straightforward dismissal of a role of 
Anderson localisation in the magnetoresistance of the manganites. Moreover, the value Vi « 2.1t 
only takes into account the random distribution of potential sources, namely RE^^ and AE"^^ ions. 
The presence of RE/AE ionic size mismatch is expected to enhance considerably the disorder effects 
|42l444 | . Note also that, in this model the critical value of disorder does not vary monotonically 
with X and shows a maximum at around x ^ 0.1. One should bear in mind that, in this model, 
changing x also changes the distribution of site energies [equation. and so Vi does not, by 
itself, characterise the disorder. 

A final remark regarding the relevance of the model used in this work is in place. The 
model given by equation ([2]) with correlated on-site disorder as in ([5]) is certainly incomplete: it 
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neglects electron-lattice coupling, orbital degrees of freedom, anti-ferromagnetic exchange between 
localised t^g spins, and electron-electron interactions between Cg electrons. Consequently, it does not 
distinguish between the three main manganite groups (large-, intermediate-, and low-bandwidth) 
and cannot reproduce many of the complex phases appearing in each of their phase diagrams 
01 • However, it is general enough and not tremendously complex to address the question of how 
important is substitutional disorder in manganites; this is the motivation for using it. 



4. Conclusions 



We have shown that a careful analysis of the Mn-site energies arising from random distribution 
of and RE^^ ions in manganites produce a probability distribution with a RMS deviation 

a « 4.6t. This RMS is a consequence of a parametrisation of screened Coulomb energies for 
which an energy shift Vi « 2. It is assumed to show up in a Mn-site whenever a NN replacement 
RE^~^ — > AE'^'^ occurs. Such a RMS value already places the system in the disorder window for 
which a MIT occurs when the ferro-paramagnetic transition takes place. We have developed a DE 
model with realistic parametrisation of on-site disorder which shows to undergo a MIT at the ferro- 
paramagnetic transition for Vi ~ M. This Vi value is slightly larger than the expected Vi ~ 2.1t, 
but it is sufficiently close to unveil the importance of substitutional disorder in manganites and show 
that disorder must be considered at least on the same foot as the coupling to the lattice. These 
findings give support to theories where disorder is a key ingredient [32], as the phase-c omp etition 
scenario for CMR, where diagonal disorder has been found to play a crucial role [Big. Such 



a key role played by disorder was also recently observed in a series of experiments in half-doped 
manganites [i^,!!^, where A-site ordered and disordered systems have been successfully prepared 
and compared. 



Acknowledgments 

The authors acknowledge financial support from POCI 2010 via project PTDC/FIS/64404/2006. 
EVC additionally acknowledges the Juan de la Cierva Program (MCI, Spain). 

References 

[1] Coey J M D, Viret M and Von Molnar S 1999 Adv. Phys. 48 167 - 293 

[2] Loktev V M and Pogorelov Y G 2000 Low Temp. Phys. 26 171 - 193 

[3] Salamon M B and Jaime M 2001 Rev. Mod. Phys. 73 583 - 628 

[4] Dagotto E, Hotta T and Moreo A 2001 Phys. Rep. 344 1 - 153 

[5] Edwards D M 2002 Adv. Phys. 51 1259 - 1318 

[6] Dagotto E 2005 Science 309 257 - 262 

[7] Zener C 1951 Phys. Rev. 82 403 - 405 

[8] Anderson P W and Hasegawa H 1955 Phys. Rev. 100 675 - 681 

[9] de Gennes P G 1960 Phys. Rev. 118 141 - 154 
[10] Kogan E M and Auslender M I 1988 phys. stat. sol. (b) 147 613 - 620 
]11] Auslender M I, Kogan E M and Tretyakov S V 1988 phys. stat. sol. (h) 148 289 - 295 



i REFERENCES 

12] Varma C M 1996 Phys. Rev. B 54 7328 - 7333 

13] Kogan E, Auslender M and Kaveh M 1999 Eur. Phys. J. B 9 373-376 

14] Li Q, Zang J, Bishop A R and Soukoulis C M 1997 Phys. Rev. B 56 4541 

15] Sheng L, Xing D Y, Sheng D N and Ting C S 1997 Phys. Rev. B 56 R7053 - R7056 

16] Sheng L, Xing D Y, Sheng D N and Ting C S 1997 Phys. Rev. Lett. 79 1710 - 1713 

17] Zarand G, Moca C P and Janko B 2005 Phys. Rev. Lett. 94 247202 

18] Pickett W E and Singh D J 1997 Phys. Rev. B 55 R8642 - R8645 

19] MiUis A J, Littlewood P B and Shraiman B I 1995 Phys. Rev. Lett. 74 5144 - 5147 

20] MilHs A J, Sliraiman B I and Mueller R 1996 Phys. Rev. Lett. 77 175 - 178 

21] Roder H, Zang J and Bishop A R 1996 Phys. Rev. Lett. 76 1356 - 1359 

22] Bozin E S, Schmidt M, DeConinck A J, Paglia G, Mitchell J F, Chatterji T, Radaelli P G, 
Proffen T and BilUnge S J L 2007 Phys. Rev. Lett. 98 137203 

23] Moreo A, Mayr M, Feiguin A, Yunoki S and Dagotto E 2000 Phys. Rev. Lett. 84 5568 - 5571 

24] Burgy J, Mayr M, Martin-Mayor V, Moreo A and Dagotto E 2001 Phys. Rev. Lett. 87 277202 

25] Burgy J, Moreo A and Dagotto E 2004 Phys. Rev. Lett. 92 097202 

26] Motome Y, Furukawa N and Nagaosa N 2003 Phys. Rev. Lett. 91 167204 

27] Motome Y, Furukawa N and Nagaosa N 2005 Lecture Notes in Physics vol 678 ed Donath M 
and Nolting W (Springer- Verlag) p 71 

28] Kumar S and Majumdar P 2006 Phys. Rev. Lett. 96 016602 

29] Sen C, Alvarez G and Dagotto E 2007 Phys. Rev. Lett. 98 127202 

30] Sen C, Alvarez G, Aliaga H and Dagotto E 2006 Phys. Rev. B 73 224441 

31] Pickett W E and Singh D J 1996 Phys. Rev. B 53 1146 - 1160 

32] Narimanov E E and Varma C M 2001 Phys. Rev. B 65 024429 

33] Mackinnon A and Kramer B 1981 Phys. Rev. Lett. 47 1546 - 1549 

34] Mackinnon A and Kramer B 1983 Z. Phys. B - Condens. Matter 53 1 - 13 

35] Kramer B and MacKinnon A 1993 Rep. Prog. Phys. 56 1469-1564 

36] Coey J M D, Viret M, Ranno L and Ounadjela K 1995 Phys. Rev. Lett. 75 3910 - 3913 

37] Boris A V, Kovaleva N N, Bazhenov A V, van Bentum P J M, Rasing T, Cheong S W, Samoilov 
A V and Yeh N C 1999 Phys. Rev. B 59 R697 - R700 

38] Salafranca J and Brey L 2006 Phys. Rev. B 73 214404 

39] Haydock R 1980 Solid State Physics vol 35 ed Ehrenreich H, Seitz F and TurnbuU D (New 
York: Academic Press) p 215 

40] Slevin K and Ohtsuki T 1999 Phys. Rev. Lett. 82 382 - 385 

41] Xiong S J, Xing D Y, Evangelou S N and Sheng D N 2003 Phys. Lett. A 311 426 - 431 

42] Sheng L, Sheng D N and Ting C S 1999 Phys. Rev. B 59 13550 - 13553 

43] Alonso J L, Fernandez L A, Guinea F, Laliena V and Martin Mayor V 2002 Phys. Rev. B 66 
104430 

[44] Kumar S and Kampf A P 2008 Phys. Rev. Lett. 100 076406 



REFERENCES 9 

[45] Akahoshi D, Uchida M, Tomioka Y, Arima T, Matsui Y and Tokura Y 2003 Phys. Rev. Lett. 
90 177203 

[46] Mathieu R, Akahoshi D, Asamitsu A, Tomioka Y and Tokura Y 2004 Phys. Rev. Lett. 93 
227202 



